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ABSTRACT 


Likelihood  ratio  statistics  for  (i)  testing  the  homoge¬ 
neity  of  a  collection  of  multinomial  parameters  against  the 
alternative  which  accounts  for  the  restriction  that  those  para¬ 
meters  are  starshaped  (cf.  Shaked,  Ann.  Statist.  (1979)),  and 
for  (ii)  testing  the  null  hypothesis  that  this  parameter  vector 
is  starshaped,  are  considered.  For  both  tests  the  asymptotic 
distribution  of  the  test  statistic  under  the  null  hypothesis  is 
a  version  of  the  chi-bar-square  distribution.  Analogous  tests 
on  a  collection  of  Poisson  means  are  also  found  to  have  asymp¬ 
totic  chi-bar-square  distributions. 

AMS  1970  Subject  Classifications:  Primary  62F05 

Secondary  62E15 

Keywords  and  phrases:  starshaped  orderings,  multinomial 

Poisson,  likelihood  ratio  tests,  chi-bar-square  distri¬ 
butions  ,  maximum  likelihood 
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1.  Introduction  and  Summary. 

Shaked  (1979)  derived  the  maximum  likelihood  estimate  of 
a  vector  of  Poisson  (normal)  means  subject  to  the  restriction 
that  this  vector  is  "starshaped."  A  vector  6  =  (01 , 02 , • • • , 0^) 
is  said  to  be  lower  starshaped  provided 

e-,  +  e„  e1  +  02  +  •  •  •  +  e,_ 


ei  +  02 

ei  *  -S~-  * 


i  0  with  an  analogous 


restriction  defining  an  upper  starshaped  vector.  Starshaped 
vectors  arise  naturally  in  reliability  theory  as  well  as  in  cer 
tain  situations  where  finite  populations  are  amalgamated.  We 
refer  the  interested  reader  to  Shaked  (1979)  for  examples  of 
parameter  sets  which  might  be  known  or  suspected  to  satisfy 
such  a  restriction. 

In  Section  2  we  consider  a  sampling  situation  where  the 
result  of  each  trial  of  our  experiment  must  be  a  member  of  a 
set  of  mutually  exclusive  events  with  corresponding  probabili¬ 
ties  Pi >P2  *  *  '  *  ,pk’  The  max^-mum  likelihood  estimate  of  the 
vector,  p  =  (p,  ,p0 , • • •  ,p,  ),  subject  to  the  restriction  that 
it  be  lower  starshaped,  is  derived  (this  sampling  situation  was 
not  considered  by  Shaked  (1979))*  This  derivation  Is  quite 
direct  and  elegant  in  light  of  the  complexity  of  the  analysis 
In  Shaked  (1979)  and  in  light  of  the  difficulties  Involved  In 
the  related  problem  of  finding  the  maximum  likelihood  estimate 
of  £  subject  to  the  restriction  p^  *p2  *•••  *Pk  (cf.  Barlow 
Bartholomew,  Bremner  and  Brunk  (1972)).  In  fact,  proofs  that 
various  algorithms  for  the  solution  to  the  latter  problem 


yield  the  desired  result  are  usually  by  induction. 


In  addition,  asymptotic  distribution  theory  for  the  like¬ 
lihood  ratio  test  of  the  homogeneity  of  Pi>p2,‘'‘>pk  aSalnst 
the  alternative  that  £  is  starshaped  and  for  testing  that  p  is 
is  starshaped  as  a  null  hypothesis  is  also  presented  in  Section  2. 
Again  the  derivations  are  relatively  direct.  In  both  situa¬ 
tions,  the  tail  probabilities  under  the  null  hypothesis  of  this 
asymptotic  distribution  turn  out  to  be  of  the  form 


**-!<*>  *  (I)1"1 


2 

where  denotes  a  standard  chi-square  variable  with  i 

degrees  of  freedom.  A  somewhat  similar  distribution  is  en¬ 
countered  in  the  problem  of  testing  homogeneity  when  the  al¬ 
ternative  is  restricted  by  Pi 51  p2  *  "  '  '  ^pk  Ghacko 

(1966))  and  for  testing  p.^  ap^  ^  ^Pk  as  a  null  hypothe¬ 
sis  (cf.  Robertson  (1978)  and  related  results  in  Robertson 
and  Wegman  (1978)).  Such  weighted  chi-square  distributions 
are  encountered  in  many  order  restricted  inference  problems 
(cf.  Barlow  et  al.  (1972)).  They  were  first  encountered  by 

Bartholomew  (1959)  and  are  usually  called  chi-bar-square  dis- 
_2 

tributions  (x  )• 

In  Section  3  we  assume  independent  samples  from  each  of  k 
Poisson  populations.  The  analysis  in  Section  2  together  with 
the  well  known  fact  that  the  joint  distribution  of  independent 
Poisson  random  variables  conditioned  on  the  value  of  the  sum 


V 


multinomial,  is  used  to  derive  maximum  likelihood  estimates 
under  the  starshaped  restriction  on  the  parameter  values. 
This  derivation  is  substantially  easier  than  the  original 
derivation  of  these  estimates  by  Shaked  (1979).  Asymptotic 
distribution  theory  for  likelihood  ratio  statistics  used  for 
testing  homogeneity  versus  starshaped  and  for  testing  star¬ 
shaped  as  a  null  hypothesis  is  also  presented. 


2.  Multinomial  Problem. 

Suppose  we  have  n  independent  trials  of  an  experiment, 

the  outcome  of  which  must  be  one  of  k  mutually  exclusive 

events  with  corresponding  probabilities,  p^  >Pp  *  "  *  '  ,r>k 

(1^  p.  *1).  Our  first  task  is  to  find  the  maximum  likelihood 
i-l‘  1 

estimate  of  the  vector  £  under  the  restriction: 


pl+P2 

(2.1)  H1  :  px  *  x2  ■  * 


P1+P2+' 


•+p 


k-1 


k-1 


We  define  a  one  to  one  transformation  of  the  parameter  space 
by  introducing  new  parameters  ®i > ®2 ’ ‘ ’ 0k-l  where 


(2.2)  6,  *  (I*  d.)/(I^p.)  :  1  =  1 ,2 ,  •  •  •  ,k-l 

1  J*1  J  j=l  J 

<P1  ■  P1  "  (1-ei-l)nj^V  1  -2. 3.  •".k-1.  Pk  *  U-Vl 


In  terms  of  the  0’s  the  likelihood  function  can  be  written 
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y 

i 

5 

J 

r 

►' 

I 

f 
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v-1  nE*  p .  np. - 

(2.3)  L(e)  =  IT_1  e1  J_1  J(i-ei)  1  ,  o^e^i, 

where  p^  is  the  relative  frequency  of  the  event  having  prob¬ 
ability  :  i  ■  1,2, • • * ,k.  The  restriction  (2.1)  becomes 


:  61  *  i/(i+l)  :  i  *1,2,**  •  ,k-l. 

a  b 

It  is  easy  to  find  the  maximum  of  the  function  0  (1-0) 
subject  to  0  ^  c  (0  sfi  s:l).  This  maximum  is  attained  at 
0  *  a/(a+b)  v  c,  where  V  denotes  the  larger  of  the  two  num¬ 
bers.  It  follows  that  the  maximum  likelihood  estimates  which 
satisfy  are  given  by 


(2.4)  ■§!  =  @i  v  (i/i+1)  :  i  *1,2,-  •  •  ,k-l 
where  =  (E^pj  )/(lJ**Pj )  and 

(2.5)  Pi  =  (l-ei_1)  0j  ;  i  = 1,2, • • • ,k-l 

with  0Q  *  0  and  *  l-G^^.  Restating  (2.5),  we  have 
established  the  following  theorem. 


Theorem  2.1.  The  maximum  likelihood  estimates  of  p1,p2,'**fpk 
subject  to  the  lower  starshaped  restriction  are  given  by 


5 


where  denotes  the  relative  frequency  of  the  event  having 

probability  p^ 


Turning  to  the  testing  problem,  let 


(2.6) 


Ho  :  P1  =P2  *  ‘  “  =pk  =  1/k’ 


We  let  A  denote  the  likelihood  ratio  test  statistic  for 
testing  Hq  against  -  HQ  (i.e.,  but  not  HQ)  and 

let  =  -2  In  Aq^.  It  is  convenient  to  write  TQ1  in  terms 

of  the  0's  as  follows: 


(2.7) 


=  2^=^nEj  =  iP^[ln  0i  -ln(i/i+l)] 


+  np1+1[ln(l~G1 )  - ln(l/i+l)] >. 


Using  Taylor's  Theorem  with  a  second  degree  remainder  term,  we 
expand  In  0^  and  ln(i/i+l)  about  ,  and  expand 
lnd-l^)  and  ln(l/i+l)  about  1-0^.  The  linear  terms  drop 


out  and  we  obtain 


6 


(2.8) 


■01 


2^ 

1=1 


nj:j  =  lP1  _  -  .2 


nl;j_  p  „  2 

+-JfLj-(e1  -jij) 
2a±  11  20^  1  1+i 


npi+1  -  _  2  np 

— ^  (e1-ei)2+  — 


2v 


1+1  (Z  In2 
2  (0J  ~1  +  1) 


2  Y 


1  "  '1 

where  Is  between  ^  and  §^;  0.^  is  between  0^  and 

i/(i+l);  v1  Is  between  (1-1^)  and  (1-6  )  and  y  is 
between  (1-?)^)  and  l/(i+l).  The  law  of  large  numbers  implies 
that,  under  HQ,  0^  converges  to  i/i+1. 

To  obtain  some  insight  into  the  asymptotic  power  of  the 
likelihood  ratio  test,  we  consider  a  sequence  of  alternatives 


satisfying 

which  converges  to  g_ 

where  p^  >  0  for 

all 

i .  We  let  £ 
~n 

denote  a  random  vector 

corresponding  to 

En’ 

i.e.  ,  np  is 

~n 

multinomial  (n,£n). 

In 

our  alternative 

parameterization  we 

have 

i 

i+1 

i 

i  +  1 

Jn,i  *  J/n.J  —  6i  ’  fV  \  PJ 


and 


i  i+1 

0  .  =  E  p  ,  /  t.  p  , 

ft > i  j=l  j=i 


Somewhat  surprisingly,  it  can  be  shown  by  conditioning  on 

i  +  1  -  -  % 

V  p_ .  ,  that  E(0  .)  =  0  •• 

L  *"3  n  ’ 1  0,1 

Since  ngn  is  multinomial  (n,pn>,  it  can  be  shown  that 


i: 


/n(pn  -  p^)  - ►  MVN(0,£)  where  £  =  is  defined  by 

pi(1  "  pi) »  i  =  j 
“PjPj  >  i  *  j  • 

(One  way  of  showing  this  is  to  verify  that  the  moments  of 

linear  combinations  of  /n(p  •  -  p  •)  converge  to  the 

n , i  n,i 

ments  of  linear  combinations  of  an  appropriate  MVN  vector,  and 
then  employ  Theorem  B  of  Serfling  (1980)  and  the  Cramer-Wold 
device . ) 

.  .  .  k  k-1 

Moreover,  if  we  define  the  function  g  =  (g  ,  •••jg^  ^):  R  — ►R 


mo- 


by 


1  +  1 


g  (x)  =  Lx./  L  X 
1  1  J  1  J 


Theorem  4.2  of  Kepner  (1979)  implies  that 

-g(£n)l  “  MVN(0,DZD') 


where 


8 


Careful  calculation  reveals  that  DID'  =  ( )  is  given  by 


(2.9) 


pi+i(^pj) 

/i+l  \ 3  ’ 

(  l  pJ 


>  i  /  J  , 


so  that,  fortunately,  we  have  asymptotic  independence  among  the 
Sn>1’s. 

If  we  recall  that  0  ,  =  0  .  v  we  can  express  the 

■Tl  )  X  24  y  1  1  *  J- 

likelihood  ratio  test  statistic  as 


(2.10)  =  I  (X  ,  +  6  ,  )2a  ,-(X  .  +  6  .)\  ,Irv  .. 

01  1=1V  n,i  n,i'  n,i  n,i  n,i'  n,i  [  Xn  A+  6n  ^sO  ] 


where 


(2.11) 


X  .  =  */n(§  -0  ,  )  ^ i+  - 3- 

n,i  n,i  n ,i  ik 


6  ,  -  ^(0  ,  -4r) 

n,i  v  n,i  i+l  ik 


i  p 

I  -4 

j-i  K 


nj.+ll  r (i  +  l)3 


y  Pn,l+ll  [W1  1 

u  2  2 

J=1  n,i  vn ,i  .  L  ik 


«s  0  for  a]  1 


Since  p  satisfies  H,  for  all  n,  6  . 

*n  1  n ,  1 


n,i. 

If 

®n,  i 

-  6i 

(finite)  as  n 

— ►  ®,  then 

6n,i 

- -4 

ITT’  so 

that 

6  « 
n,i 

p  l 

i+1* 

Recalling  how  a  . , 
n ,  i 

Bn.i 

»  Yn 

,  ,  and 

Vn,i 

were 

obtained , 

it  then  follows  that 

a  , 

P  , 

1  and 

b  ,  - 

-P-  1 

as  n  — * 

n,i 

n,i 

In  this  situation,  X  ,  - *  Z,  ,  where  Z  is  a  n(0,l) 

n  ,1  1  i 

random  variable.  Using  Theorem  4.9  of  Billingsley  (1968),  we 
have 


ft 

(X  ,  +6  ,a  .  ,b  ,)  - ►  (Z,  +6,, 1,1). 

n,i  n,i*  n,i’  n,i  i  i*  * 


■?  1 

Then  noting  that  the  function  h  :  RJ  — *  R  defined  by 


2  2 

h(x,y,z)  =  x  y  -x  z  Ij-xs0] 


is  continuous,  we  may  use  Theorem  5.1  of  Billingsley  to  say 


(X  .  +  6  ,)2a  ,  -  (X  ,  +  6  ,)2b  ,  Irv  .. 

n,I  n,i  n,i  n,i  n,i  n,i  [X  , +6  , sOJ 

n  j  i  n  j  _L 


3 


(Zi+6i}  -  (W  I[Z1  +  61*0] 


-  C(Zi+61)  VO]' 


In  the  event  that  6  ,  — ►  ®,  It  can  be  shown  that  a  ,  is 

n,I  n  ,1 

bounded  away  from  zero  asymptotically  while  X  ,  converges 

n ,  i 


10 


in  distribution.  Thus 


(Xn,i  +  6n,i)  an,i  "  (Xn,i  +  6n,i)  bn,i  I[X  ,+6  *0]  P_>  “• 

n  j  i  n  j  i 

We  have  thus  established  the  following  theorem. 

Theorem  2.2.  If  satisfying  converges  to  g(p^>0  for  all  i)  , 

and  if  6n  ±  (as  defined  in  (2.11))  - *  61  (possibly  •)  for 

i“l»***,k-l,  then  Tq  ^  is  distributed  asymptotically  as 


(2.12) 


U  =  Z  [  ( Z .  +  6.  )  v  0]  , 
i=l  1  1 


where  are  independent  n(0,l)  random  variables. 

Of  course  the  distribution  of  the  random  quantity  In 
(2.12)  is  very  intractable,  except  under  the  null  distribution 
Hq  (*±  =0,  i=l,***,k-l)  when  it  becomes  surprisingly  nice. 

To  elaborate,  suppose  I  is  a  subset  of  {l ,2 , • • •  ,k-l }  and 
let  Ej  be  the  event  Ej  =  [Z1  2=0;  i  €  I  and  Z±  <0;  i  £1]. 
Then,  for  any  real  number  u. 


P[U  >u,Ej]  =  P[Z1€I  Z^au,  Z^O;  ifl,  Z±  <0;  i  JEI] 

=  P[Z1€I  z\  2u,  Z1  20;  i  €I]*P[Z1  <0;  1  jf  H 
*  P[Z1€1  Z^  2u  |Z1  20;  i  €  I]* (l/2)k-1 


P[*£  *u]-(l/2)k_1 


11 


where  m  is  the  number  of  elements  in  I.  The  last  step  fol¬ 
lows  from  Lemma  B  on  page  128  of  Barlow,  Bartholomew,  Bremner 
and  Brunk  (1972).  Partitioning  the  event  [Biu]  by  inter¬ 
secting  it  with  all  such  events,  Ej,  we  obtain  the  expression 
for  P[U  2s  u]  given  in  the  following  theorem. 

Theorem  2.3-  IT  HQ  is  true  then 

Umn_P[T01^]  =  ^(kl1)(l/2)k_1  PCX^t]  - 

p 

for  all  real  t  (x0E0). 

The  expression  (2.12)  is  useful  in  getting  a  feeling  for 
the  asymptotic  power  of  the  test  compared  to  the  usual  unre¬ 
stricted  test  for  homogeneity  of  multinomial  parameters.  In 
particular,  under  the  conditions  of  Theorem  2.2,  the  asymptotic 
distribution  of  -2  In  A  where  A  is  the  unrestricted  likeli¬ 
hood  ratio  (or  of  the  usual  Pearson  chi-square  goodness  of  fit 
test)  is  the  same  as 

k-1  0 

U'  -  L  (Z.  +  6,  r 
1  1  1 

where  ^i»*‘'>Zk-l  are  *nc*ePendent  n(0,l)  random  variables. 
Under  Hq,  U'  is  x  (k-1)  and  hence  by  looking  at  Theorem 
2.3,  we  see  that  its  critical  point  must  be  substantially 


larger  than  for  testing  Hq  against  H^.  However,  as  the 
6^ ’ s  become  larger,  more  nonzero  terms  enter  into  (2.12),  so 
that  U  and  U'  become  more  nearly  equivalent.  The  smaller 
critical  point  of  the  restricted  test  implies  that  its  power 
must  be  larger  eventually  than  that  of  the  unrestricted  test. 

We  now  turn  to  the  problem  of  testing  as  a  null  hypoth¬ 
esis.  Since  the  unrestricted  maximum  likelihood  estimate  of 
0^  is  equal  to  0^,  it  follows  directly  by  writing  the 
likelihood  ratio  in  terms  of  0  and  ¥  and  expanding  In  ¥1 

(lnd-0^)  about  0i  ((1-6^))  that  our  test  statistic  can  be 

r  i 


written  as 


-2  In 


i 


n  ’  ^i_6i^ 


where  is  between  0i  and  3^  (and  thus  converges  a.s.  to 

_  * 

0^)  and  is  between  1-0^  and  1-6^  (and  thus  converges 

a.s.  to  1-0^. 

By  employing  arguments  similar  to  those  used  in  Theorem 
2.2,  we  are  led  to  the  following  theorem.  (Note  that  we  do  not 
need  to  restrict  to  H^.) 


Theorem  2.H.  If  p  converges  to  p  and  5  .  (as  defined 

~n  ~  11,1 

in  (2 . 11))  converges  to  6^  ( are  possible  values)  for 

i  *  1, * ' * ,k-l,  then  T^  is  distributed  asymptotically  as 

k-i  , 

V  =  I  [(Z.+M  AO]^ 

1=1  1  1 


where  Z^, 


k-1 


•  *  « 


are  independent  n(0,l)  random  variables. 


We  note  that  if  i-1  Ep.  >  (i+l)~'L  E  P.j  then  6  . 

^  J  1  J  II  ,  i 

t  h 

In  this  case,  the  ib  term  in  (2.13)  is  zero  and  can  be 
ignored,  leading  to  the  following  theorem. 


-1 


i+1 


Theorem  2.5.  If  p  satisfies  H1  then 


-2. 


p['r1  =  ^(t) 


where  m  is  the  number  of  subscripts,  i,  such  that 
i  =  (i+1)  ^^j-^Pj’  i  = 1,2,  • • • ,k-l.  In  addition. 


llmn-*“  Pp[TlIt:l  *  llnW  PH0tTl 

where  P„  [Tn  it]  is  the  probability  of  the  event  [T,  it] 

Hq  i  1 

computed  under  HQ • 

We  note  that  Theorems  2.2  and  2.4  imply  that  the  likeli¬ 
hood  ratio  tests  considered  here  are  consistent  in  the  sense 
that  for  p  lying  in  the  region  defined  by  the  alternative 
hypothesis,  the  power  function  must  converge  to  one. 


3.  Poisson  Problem. 


Suppose  we  have  a  random  sample  of  size  n  from  each  of 
k  Poisson  populations  having  means  *•]_  >  ^2  ’  ’  "  ’  5  Xk '  Shaked 
(1979)  found  the  maximum  likelihood  estimate  of 
X=  (X^,X2, • • • ,Xk)  subject  to  the  restriction  requiring 
that  X  is  lower  starshaped: 


(3.1) 


X  +X 

X  2  X  <  > 
*12 


Xl+x2+"  ’  +  Xk 
2  — - - -  2  0. 


This  result  can  be  found  in  a  straightforward  fashion  using  the 
results  in  Section  2  together  with  the  fact  that  the  conditional 
distribution  of  independent  Poisson  variables,  given  their  sum, 
is  multinomially  distributed. 

We  first  write  the  likelihood  function  in  terms 
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of  the  variables  » ^2  *  *  *  * ’ where 

k  k 

(3-2)  j =  X^/L^  X^;  i  =  1 ,2  ,  •  ■  *  ,k-l  and  #k  = 

(Xi  =  **i0k;  1  =  !>2>  *  * '  »k-l»  ^  =  ^k  ~  ^i=i0iiZik)  •  The  restric“ 

tion  that  X  is  starshaped  is  equivalent  to  requiring  that 

^  is  starshaped,  or 


0-1  +0p 

(3.3)  H * 


(k-l)"1I^“^01  *  k 


-1 


and  these  restrictions  do  not  involve  0^.  The  likelihood  func¬ 
tion  is  proportional  to 


(3-4) 


-n0. 


£h=ixn 


_  t  h 

where  is  the  mean  of  the  sample  from  the  i  population; 

ial,2,**',k.  Because  does  not  restrict  0^,  the  two 

factors  in  brackets  may  be  maximized  independently.  Using  the 

results  from  Section  2  on  the  first  factor  and  an  easy  analysis 

involving  the  derivative  of  the  second  factor,  we  obtain  the 

restricted  maximum  likelihood  estimates  as  follows: 


(3.5) 


?i  =  ?i  V  (i/i+l)  :  i  =  1,2,  •  •  •  ,k-l 


0^.  (Note  that 


where  0^  *  V^-1  xj  •  and  -k  *  $.1  xj 
^1*^2*“  *  are  unrestr^-cted  maximum  likelihood  estimates 
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of  *  *  *  *  ^  Using  the  Invariance  property  of  maximum 

likelihood  estimation,  we  have  the  following  theorem. 

Theorem  3»1  (Shakgd).  The  maximum  likelihood  estimates  of 
X 1 , X2 , • • • , Xk  subject  to  the  restriction  are  given  by 


(3.6)  and 


Xi  =  £*  ±-  V  ill  ‘Lj=iXj;  1=1,2 

L>ixj  J  J 

i 

(I^-A)[1-C(?\viir) 

L  >n  J 


•,  i  =1,2,-  •  •  ,k-l 


The  likelihood  ratio  statistic  for  testing  HQ  : 

X^  *  X2  * • • •  =  X^  against  the  alternative  H^-HQ  can  be  written 


in  terms  of  the  jzs’s,  as  follows 


(3-7) 


(l/k) 


n5i-ixi 


k-1  -nX,  (  yk-l  -  \  nXv 
“i-i  *  1  -  ri=i  *1  k 


If  we  let  S,,,  =  -2  In  and  let  Y  =  nZ  X. ,  then  given 

01  01  i = i  I 

that  Y  =  y,  the  joint  conditional  distribution  of 

y^l,y^2’' *  * ,y(1  is  multinomial  with  parameters  y  arid 

x,  /Z^  x, ,  x  /Z  x, ,  ■  •  * ,  x.  /Zk  X,  . 

1  i-1  i  2  1*1  i*  k  i*i  i 
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If  we  let  \n  satisfying  converge  to  X  ( X^  >  0 

for  all  1)  such  that 


(3.8) 


(  1  1+1  4  \ ' 

=  I  X  ./EX  ,  -txt) 

Vj=1  n,J  J=1  nj  1+1/ 


(1+1) 


6^  (possibly  ®) 


and  let  X  .  denote  the  corresponding  independent  sample 

n,i 

means  which  occur  in  Sq^  =  -2  In  /q|>  then  using  the  Domi¬ 
nated  Convergence  Theorem: 


lim  P(S^}  at) 
n-»®  ux 


=  E[P(U  st)] 


lim  ECPrs^^  st  |Y  )) 
_  01  1  n 

n-*® 


E[lim  P(S<"}  st  |Y  )] 
n  -*00  u  1  n 


=  P[U2  t] 


where  U  is  distributed  as  in  Theorem  2.2.  Thus  SQ1  is 

distributed  asymptotically  exactly  as  TQ1  in  Section  2  (with 

p  .  replaced  by  X  .  in  the  6  . 's).  In  particular,  if 

n,i  n  j l  n  j l 

Hq  is  true,  the  asymptotic  tail  probabilities  are  given  as 
in  Theorem  2.3. 
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Theorem  3.2.  If  Hq  is  true,  then 

P^-S01  =  tyc-l^' 

Let  us  now  consider  the  problem  of  testing  as  a  null 

hypothesis.  The  likelihood  ratio  A  can  be  written 


since  =  0R. 

If  we  let  S1  =  -2  In  then  the  same  type  of  reasoning 

used  in  the  previous  argument  can  be  applied  to  show  that 

(3-10)  lim  P(S  a  t)  =  lim  E(P(S,  *t  |Y)) 

n-*«  x  n-*® 

=  EClim  P(S  at t  |Y)  ] 

n  *♦* 

=  p[V  st] 

where  V  is  defined  in  Theorem  2.4  and  6  .  is  defined  as 

n,i 

in  (3.8).  It  follows  from  (3.10)  that  if  X  satisfies  II,, 

i**  JL 

then  the  itfl  term  of  V  goes  to  zero  if 
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-l  1  i1+1 

1  Z  x  >  (1+1)"1  Z  x 

3-1  J  3=1  J 


which  leads  to  the  following  theorem. 


Theorem  3-3.  If  \  satisfies  H^,  then 


lim  PXCS1  &  t]  =  X^(t) 


n-+ 


where  m  is  the  number  of  distinct  i  such  that 


-i  1  li  +  1 

i  I  X  =  (i+1)"1  Z  X 
3=1  J  3=1  J 

Moreover 


(3.11)  sup  lim  P,[S  it]  =  P„  [S.  st]  =  y?  ,  (t) 
X€H,  n-*»  ~  1  H0  1  k_1 


where  ( • )  is  computed  under  the  assumption  that  is  true. 


Note  that  Theorem  3.3  enables  us  to  construct  likelihood 


ratio  tests  of  a  particular  size  asymptotically  when  testing 
H1  versus  all  other  alternatives.  Of  course  (3.8)  and  (3-10) 
assure  us  that  our  tests  are  asymptotically  consistent  in  the 
sense  that  if  X  is  in  the  region  of  the  alternate  hypothesis, 

A/ 

the  probability  of  rejecting  the  null  hypothesis  converges  to 
one  as  n  — ♦  ®. 

It  should  be  noted  that  even  though  Shaked  (1979)  allows 
the  more  general  starshaped  ordering 


2  2  3 

Xn  a  £w, X  /£w,  a  £ w 
-L  i  J  J  i  J 


jviki 


k  k 

£w, X,/£w,  £ 
1  J  J  n  J 


o. 


his  restriction  that  the  sample  size  from  the  1th  population 
be  proportional  to  effectively  reduces  the  problem  to  the 
one  considered  earlier. 


4.  Concluding  Remarks. 

F.  T.  Wright  called  our  attention  to  the  work  of  Shaked 
(1979)  after  we  had  carried  through  much  of  the  research  in 
this  paper.  Actually  our  original  analysis  neglected  the 
nonnegativity  restriction  and  we  termed  the  restriction 
"decreasing  on  the  average."  More  specifically  we  should  have 
termed  It  "decreasing  on  the  average  from  the  left"  since 

e,  +  e5  8,  +  e5  +  e,  e,  +  e?  + 

8^  a  —=~2 — —  2  — - 1 ^  is  not  equivalent  to  - ^ - - 

6,  +  e? 

— ■  a  8 ^  (i.e.,  increasing  on  the  average  from  the  right). 


A 


a 
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It  is  clear  that  the  restrictions  increasing  on  the  average 
from  the  right,  increasing  on  the  average  from  the  left  and 
decreasing  on  the  average  from  the  right  can  be  handled  by 
analysis  similar  to  that  in  Sections  2  and  3. 

The  phrase  "decreasing  on  the  average"  also  calls  to  mind 
the  restriction 

H2  :  i-1  ^=16J  *  (k-i)-1  £j=1  +  16j  5  i  =  1,2,  •  •  •  ,k-l. 

An  equivalent  way  of  stating  is  i-1Ij=16j  a  k-1L^=10j  ; 

i  = 1,2, • • •  ,k-l.  We  note  that  the  order  restrictions  specified 
in  H2  are  less  restrictive  than  those  imposed  by  which 

in  turn  are  less  restrictive  than  0.^  £  ®i  +  i»  i=l**,'»k-l. 

In  the  multinomial  setting,  maximum  likelihood  estimates  of  p 
subject  to  H2  and  distribution  theory  for  testing  HQ  vs. 
H2-Hq  and  for  testing  H2  as  a  null  hypothesis  can  be  found 
in  Robertson  and  Wright  (1980).  Again,  the  asymptotic  distri¬ 
bution  is  a  chi-bar-square. 

The  10$,  5%  and  1%  cutoff  values  of  x2 ( ' )  for 
k»2,3»***>15  are  given  in  Table  1. 
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Table  1. 

Critical  values  for  xk-1(t)  =  P[x^t] 


a 


k  > 

.  10 

.05 

.  01 

3 

2.95 

4.23 

7.28 

4 

4.01 

5.44 

8.77 

5 

4.95 

6.50 

10.02 

6 

5.84 

7.48 

11.18 

7 

6.67 

8.41 

12.26 

8 

7.48 

9.29 

13.31 

9 

8.26 

10.15 

14 . 29 

10 

9.02 

10.99 

15.29 

11 

9.76 

11.79 

16.21 

12 

10.49 

12.59 

17-12 

13 

11.22 

13.38 

18.01 

14 

11.93 

14.15 

18.91 

15 

12.63 

14.91 

19.78 
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